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The Boltzmann equation for inelastic Maxwell models is considered to determine the rheological 
£Nj , properties in a granular binary mixture in the simple shear flow state. The transport coefficients 

(shear viscosity and viscometric functions) are exactly evaluated in terms of the coefficients of resti- 
tution, the (reduced) shear rate and the parameters of the mixture (particle masses, diameters and 
concentration). The results show that in general, for a given value of the coefficients of restitution, 
the above transport properties decrease with increasing shear rate. 
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I. INTRODUCTION 
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CJ ■ It is well-recognized that granular matter can be modeled by a fluid of hard spheres with inelastic collisions. In 
the simplest version, the grains are assumed to be smooth so that the inelasticity is only accounted for by a constant 
coefficient of normal restitution. For sufficiently low-densities, the (inelastic) Boltzmann equation has been solved by 
means of the Chapman-Enskog method [l[ and the Navier-Stokes transport coefficients have been obtained in terms 
of the coefficient of restitution 0]. Moreover, some analytical results in far from equilibrium situations have been 
also reported in the literature for inelastic hard spheres However, due to the complex mathematical structure of 
the Boltzmann collision operator, it is generally not possible to get exact results from the Boltzmann equation for 
inelastic hard spheres and consequently, most of the analytical results have been derived by using approximations 
£h and/or kinetic models. 

I ■ As in the elastic case, a possible way to partially overcome the above difficulties is to consider the so-called inelastic 
Maxwell models (IMM), where the collision rate is independent of the relative velocity of the two colliding particles B. 
q Thanks to this property, nonlinear transport properties can be exactly obtained in some particular problems 0, 

' for IMM without introducing additional and sometimes uncontrolled approximations. In addition, apart from their 
'— academic interest, it has also been shown that in some cases the results derived from IMM compare well with those 
obtained for inelastic hard spheres Q and even recent experiments Q for magnetic grains with dipolar interactions 
^ , are well described by IMM. All these results stimulate the use of this interaction model as a toy model to characterize 
the influence of the inelasticity of collisions on the physical properties of granular fluids. 
\ The aim of this paper is to determine the rheological properties (shear stress and normal stress differences) in a 
CNj ■ binary mixture described by the Boltzmann equation for IMM and subjected to the simple or uniform shear flow 
(USF). This state is perhaps one of the most widely studied states in granular gases 1. At a macroscopic level, the 
. USF is characterized by constant partial densities n r , a uniform granular temperature T, and a linear velocity profile 
7-H ' u x = ay, where a is the constant shear rate. Under these conditions, the mass and heat fluxes vanish by symmetry 
ON . and the pressure tensor Py is the only relevant flux in the problem. Conservation of momentum implies P^ y = const 
■ while the energy balance equation reads 



~T=--P X y-TC (1) 



d T _ 2a 

' dt dn 

h : 

>t i where d is the dimensionality of the system (d = 3 for spheres and d = 2 for disks) and C is the inelastic cooling rate. 
Equation fl} clearly shows that the temperature changes over time due to two competing effects: the viscous heating 
term aP xy and the inelastic collisional cooling term £T. While the first term is inherently positive (since P xy < 0), the 
second term is negative since Q > 0. Depending on the initial condition, one of the effects prevails over the other one 
so that the temperature either decreases or increases in time, until a steady state is reached for sufficiently long times. 
Given that in this steady state the (reduced) shear rate is enslaved to the coefficients of restitution, the problem is 
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inhere ntly non-Newtonian (and so, beyond the scope of the Navier-Stokes description) in highly inelastic granular 
gases 

As in a previous paper for a single gas Q, the rheological properties of the granular binary mixture are given in 
terms of a collision frequency Uq, which in principle can be freely chosen. Here we will consider two classes of IMM: 
(i) a collision frequency vq independent of temperature (Model A) and (ii) a collision frequency vo(T) monotonically 
increasing with temperature (Model B). Model A is closer to the original model of Maxwell molecules for elastic 
gases [O, [HI while Model B with vq(T) oc y/T is closer to inelastic hard spheres. The possibility of having a general 
temperature dependence of v$(T) for inelastic repulsive models has been recently introduced by Ernst and co-workers 
As will be shown later, Models A and B lead to the same results in the steady state limit. In particular, the 
reduced shear rate a* = a/vo m the steady state is a universal function a*(ctij) of the coefficients of restitution ctij and 
the parameters of the mixture. However, since a* does not change in time for Model A, a steady state does not exist 
except for the specific value a* = a*(aij). Consequently, a non- Newtonian hydrodynamic regime (where a* and 
are independent parameters) is reached in the long-time limit where the combined effect of both control parameters 
on the rheological properties can be studied analytically for Model A. This is an interesting new added value of this 
simple model. 

The plan of the paper is as follows. We first introduce in Section [II] the Boltzmann equation framework for IMM; 
collisional moments are worked out. In Section ITm we will introduce driving through a macroscopic shear, and consider 
in particular the USF problem. We shall subsequently focus on rheological properties in Section HVl where our results 
for the non-linear shear viscosity and viscometric functions will be reported. Finally, conclusions will be drawn in 
Section El 



II. THE BOLTZMANN EQUATION FOR IMM. COLLISIONAL MOMENTS 



Let us consider a binary mixture of inelastic Maxwell gases at low density. In the absence of external forces, the 
set of nonlinear Boltzmann equations for the mixture reads 

(^+v-V^/ r (r,v;t)=^;j„[v|/ r (t),/.(t)] , (2) 

^ ' s 

where / r (r, v;£) is the one-particle distribution function of species r (r = 1,2) and the Boltzmann collision operator 
J rs [v|/ r , f s ] describing the scattering of pairs of particles is 



Jrs [vi|/ r ,/J - J dv 2 J da [a-V r (r,vi,i)/ s (r,v^i) - / r (r, vi,t)/ a (r, v 2 ,t)] . (3) 

Here. 

n r = J dv/ r (v) (4) 

is the number density of species r, oj rs is an effective collision frequency (to be chosen later) for collisions of type 
r — s, fid = 27r d / 2 /r(d/2) is the total solid angle in d dimensions, and a rs < 1 refers to the constant coefficient of 
restitution for collisions between particles of species r with s. In addition, the primes on the velocities denote the 
initial values {v'^v^} that lead to {vi,V2} following a binary collision: 

v'i = vi - n S r (l + a" 1 ) (ct ■ gl 2 )CT, (5) 



V2 = v 2 + firs (1 + a,,/) {a ■ gi 2 )CT , (6) 

where gi 2 = vi — v 2 is the relative velocity of the colliding pair, <x is a unit vector directed along the centers of the 
two colliding spheres, and fi rs = m r /(m r + m s ). 

The effective collision frequencies uj rs are independent of velocity but depend on space an time through its depen- 
dence on density and temperature. Here, as in previous works (7| for monocomponent gases, we will assume that 
uj rs oc n s T q , with q > 0. The case q = (a collision frequency independent of temperature) will be referred to as 
Model A while the case q > (collision frequency monotonically increasing with temperature) will be called Model 
B. Model A is closer to the original model of Maxwell molecules for clastic gases [llOj while Model B, with q = ±, 
is closer to hard spheres. 



3 



Apart from n r , at a hydrodynamic level, the relevant quantities in a binary mixture are the flow velocity u, and 
the "granular" temperature T. They are defined in terms of moments of the distribution f r as 



pu = p r u r 



Y J dvm r vf r (v), (7) 



nT = Y n r T r = Y J dv^V 2 f r {v), 



(8) 



where p r = m r n r , n = ri\ + n,2 is the total number density, p = p\ + pi is the total mass density, and V = v — u is 
the peculiar velocity. Equations ([7|) and {8} also define the flow velocity u r and the partial temperature T r of species 
r, the latter measuring the mean kinetic energy of species r. As confirmed by computer simulations (T3 |. experiments 
[lH and kinetic theory calculations [l6j . the global granular temperature T is in general different from the partial 
temperatures T r . 

The collision operators conserve the particle number of each species and the total momentum but the total energy 
is not conserved: 

j dvJ rs [v\f r ,f s ] =0, (9) 

2 2 

Y Y m - / dvvJ ™ M/r, ^ = ' ( 10 ) 



2 2 



Yll m -\ dvV 2 J rs [v\f r J s ] = -dnTC , (11) 

r=l s=l ^ 

where £ is identified as the total "cooling rate" due to inelastic collisions among all species. At a kinetic level, it is 
also convenient to discuss energy transfer in terms of the "cooling rates" £ r for the partial temperatures T r . They are 
defined as 

Cr = Y C™ = -j^Y Y j dvm r V 2 J rs {v\f ri /,] . (12) 

s r r s 

The second equality in (fT2|) defines ( rs . The total cooling rate £ can be expressed in terms of the partial cooling rates 
as 

2 

C = T- 1 Y XrTrCr, (13) 
r=l 

where x r = n r /n is the mole fraction of species r. 

From Eqs. ([9]) to (|11|) . the macroscopic balance equations for the mass, momentum and energy of the binary mixture 
can be easily obtained. They are given by 

D t n r + n r V ■ u + V ' Jr = , (14) 
m r 

Dtu + p _1 V-P = 0, (15) 

D t T -~Y — ^ + (V ■ q + P : Vu) = -(T . (16) 
77 m r an 

r—l 

In the above equations, D t = dt + u • V is the material derivative, 



jr = m r J dvV/ r (v), (17) 
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is the mass flux for species r relative to the local flow. 



Y, / rfvm r VV/ r (v), 



(18) 



is the total pressure tensor, and 



°1 = E / dv±m r V 2 Vf r (v) 



(19) 



is the total heat flux. 

The main advantage of the Boltzmann equation for Maxwell models (both elastic and inelastic) is that the (col- 
lisional) moments of J rs [fr,fs] can be exactly evaluated in terms of the moments of f r and f s without the explicit 
knowledge of both distribution functions [ll[ ■ This property has been recently exploited [2C| to obtain the detailed 
expressions for all the second-, third- and fourth-degree collisional moments for a monodispcrse gas. In the case of a 
binary mixture, only the first- and second-degree collisional moments have been also explicitly obtained. In particular 



dvm r WJ rs [f r ,f s ] = —M sr {l + a rs ) {2p s P r - (j r j s + j s j r ) 

p s d 

2 

d , . . . " 

~ [PrPs + psPr) ~Jr is 1 



(20) 



where 



dv m r VV f r 



(21) 



p r = n r T r = trP r /d is the partial pressure of species r, and 1 is the d x d unit tensor. It must be remarked that, in 
general beyond the linear hydrodynamic regime (Navier-Stokcs order), the above property of the Boltzmann collision 
operator is not sufficient to exactly solve the hierarchy of moment equations due to the free-streaming term of the 
Boltzmann equation. Nevertheless, there exist some particular situations (such as the simple shear flow problem) for 
which the above hierarchy can be recursively solved. 

The cooling rates £ rs defined by Eq. (|12j) can be easily obtained from Eq. (|20|) as 



2u> r 



■fi sr (l + a rs ) 



Msr , .0 r + 9 s fi sr (l + a rs ) - 1 . . 

2 U s dp s pr 



where 



Tn r ^ — ^ _ i 
7r V 



(22) 



(23) 



and 7 r = T r /T . Equation (|22p provides the relationship between the collision frequencies u> rs and the cooling rates 
Cr S . This relationship can be used to fix the explicit forms of u> rs . As in previous works on inelastic Maxwell mixtures 

BUI 

chosen here to guarantee that the cooling rate for IMM be the same as that of inelastic hard spheres 
(evaluated at the local equilibrium approximation) [l6| . With this choice, one gets 



012 



d-1 



1/2 



"a, = A{q)nT q , 



(24) 



where the value of the quantity A(q) is irrelevant for our purposes. Upon deriving (|24[) use has been made of the fact 
that the mass flux j r vanishes in the local equilibrium approximation. In the remainder of this paper, we will take the 
choice (|2~4"|) for lo ts . The results for IMM 0] obtained with the latter choice in the steady shear flow problem compare 
very well with those theoretically obtained for inelastic hard spheres in the first Soninc approximation and by means 
of Monte Carlo simulations [23| . 



5 



III. UNIFORM SHEAR FLOW 



Let us assume that the binary mixture is under USF. The USF state is macroscopically defined by constant densities 
n r , a spatially uniform temperature T(t) and a linear velocity profile u(y) = Ui(y) = u 2 (y) = ayx., where a is the 
constant shear rate. Since n r and T are uniform, then j r = q = 0, and the transport of momentum (measured 
by the pressure tensor) is the relevant phenomenon. At a microscopic level, the USF is characterized by a velocity 
distribution function that becomes uniform in the local Lagrangian frame, i.e., / r (r,v;£) = f r (V,t). In this frame, 
the Boltzmann equation ((2]) reads [13] 



§- t h - aV v-^rh = Jn[h,h] + MhJ 



and a similar equation for f 2 . Equation (|25[) is invariant under the transformations (V x , V y ) — > {—V x , — V y ), Vi — > —Vi, 
with i ^ x,y. This implies that if the initial state is compatible with the latter symmetry properties, then the solution 
to (|25[) has the same properties at any time t > 0. Note that the properties of uniform temperature and constant 
densities and shear rate are enforced in computer simulations by applying the Lees- Edwards boundary conditions [l7| , 
regardless of the particular interaction model considered. In the case of boundary conditions representing realistic 
plates in relative motion, the corresponding non-equilibrium state is the so-called Couette flow, where densities, 
temperature and shear rate are no longer uniform [la . [l9| . 

As said before, the rhcological properties of the mixture are obtained from the pressure tensor P = Pi + P2, where 
the partial pressure tensors P r (r — 1,2) are defined by Eq. (|2"Tj) . The elements of these tensors can be obtained by 
multiplying the Boltzmann equation (|2l)|) by m r VV and integrating over V. The result is 

d 

— P Mj + '/,/. /'1 a , + a jk P 1M + B u P 14j + B 12 P 2 ,ij = (A n pi + A 12 p 2 ) 5ij, (26) 
where use has been made of Eq. (|2H|) (with j r = 0). In Eq. (|2"6"1) . = adi X Sj y and we have introduced the coefficients 

^^(iw + ^aw, (27) 

A - = jf-27/^ 1 + a ^ (28) 

Bn = ^ (1 + an)(d + 1 - a n ) + ^ ^(1 + a 12 ) [d + 2 - ^(1 + a 12 )] , (29) 



-B12 = —3^12- 
a 



(30) 



A similar equation can be obtained for P 2 , by adequate change of indices 1 <-> 2. The balance equation |T]) for the 
temperature can be easily obtained from Eq. (|26p . In reduced units, Eq. (JTJ) can be written as 



-c 



d xyy 



(31) 



where C* = C/^Oj a * = a /f0i P£ y — Pxy/Pi P = n ^ being the hydrostatic pressure. The expression for can be 
obtained from Eqs. (|12|) . (fT3|) and (|22|) when one takes j r = 0. It is given by 



EE 

r— 1 s— 1 



37 y .Z^ 5 



0~ rs 
0~12 



d-1 



6,, + 6< 



1/2 



7 r /i sr (l + a sr ) 



-(1 + a rs )- 



+ 



(32) 



As said in the Introduction, Eq. (|3"Tj) shows that the temperature changes in time due the competition of two 
opposite mechanisms: on the one hand, viscous heating (shearing work) and, on the other hand, energy dissipation 
in collisions. Moreover, the reduced shear rate a* is the non-equilibrium relevant parameter of the USF problem since 
it measures the distance of the system from the homogeneous cooling state. It is apparent that, except for Model 
A (q = 0), the collision frequency vo(T) oc T q is an increasing function of temperature (provided q > 0), and so 
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a*(t) oc T(t)~ q is a function of time. Consequently, for q ^ 0, after a transient regime a steady state is achieved in the 
long time limit when both viscous heating and collisional cooling cancel each other and the mixture autonomously 
seeks the temperature at which the above balance occurs. In this steady state, the reduced shear rate and the 
coefficients of restitution are not independent parameters since they are related through the steady state condition 

a sPs,xy = — 2^ ' 

where we have called a* and P% xy the steady-state values of the (reduced) shear rate and the pressure tensor. On the 
other hand, when q = 0, dta* = and the reduced shear rate remains in its initial value regardless of the values of 
the coefficients of restitution a rs . As a consequence, in the case of Model A, there is no steady state (unless a* takes 
the specific value a* given by the condition (|33p) and a* and a rs are independent parameters in the USF problem. 
The analytical study of the combined effect of both control parameters on the pressure tensor is the main goal of this 
paper. 



IV. RHEOLOGICAL PROPERTIES 



In order to characterize the nonlinear response of the system to the action of strong shearing, it is convenient to 
introduce the non-linear shear viscosity rj* and the viscometric functions and "I^ as 

* I *\ 1/0 X V (OA\ 

r) ( a ) = , (34) 

p a 

u 2 P — P u 2 P — P 

y*Ja*) = ^ xx 9 vv , %{a*) = ^ zz VV . (35) 
pa pa 

The viscosity function rf (a* ) is a measure of the breakdown of the linear relationship between the shear stress P xy 
and the shear rate (Newton's law), while the first and second viscometric functions f J 2 ( a *) represent the normal 
stress differences. The explicit form of the above functions depends on the interaction model considered. 



A. Model A 



In Model A (q = 0), the collision frequency is independent of temperature and the reduced shear rate a* is a 
constant. Thus, Eq. (j26|) and its counterpart for P2 constitute a linear homogeneous set of coupled differential 
equations. In fact, it is easy to see that the relevant elements of the partial pressure tensors are the xy-elements 
along with the diagonal ones. As expected, the remaining elements tend to zero in the long-time limit. Moreover, 



from Eq. (J^BJ) is also easy to prove that, for long times, P, 



p — 

J r.zz — 



= P 



r,dd< 



Thus, according to Eq. (f3"S*)) . the 



second viscometric function ^ = 0. This is a particular property of IMM since f J ^ f° r inelastic hard spheres 



|22j, although its magnitude is always much smaller than that of As a consequence, the relevant elements of the 
partial pressure tensors are P r , X x — Pr — {d — l)P r ^ vy , P r ,yy, and P r ^ xy with r = 1, 2. 

As in the monocomponent granular case Q, one can check that, after a certain kinetic regime lasting a few collision 
times, the scaled pressure tensors P*+j = Pr,ij/p reach well-defined stationary values (non-Newtonian hydrodynamic 
regime) , which are non- linear functions of the (reduced) shear rate a* = (x/vq and the coefficients of restitution. In 
terms of these scaled variables and by using matrix notation, Eq. (|26[) can be rewritten as 



CV= Q, 

where V is the column matrix defined by the set 

■p _ r p* p* p* p* p* p* \ 
1 — l-'ljXX' r l,yy r l,xy' r 2, in r 2, yy 2,xyii 

Q is the column matrix 



(36) 
(37) 



Q 



A* n Pl + A* 12 p* 2 


A22P2 + A21P* 
A22P2 "f" A21P1 





(38) 



and C is the square matrix 
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£ = 



Here, p* = p r /p = x r j r , A* s = A rs /vn and B* 
for long times the temperature behaves as 



/ #1*1 + A 





2a* 


#12 








\ 





#n + A 








#12 








a* 


B*i + A 








#12 




#21 








i?22 + A 





2a* 







#21 








S| 2 + A 







V o 





#21 





a 


#22 + 





(39) 

= B ts Jvq. Moreover, on physical grounds it has been assumed that 
T{t) = T(0)e Al/ °* (40) 



where A is also a nonlinear function of a* , a rs and the parameters of the mixture. The (reduced) total pressure tensor 
P*j = Pij/p of the mixture is defined as 



P* 



P 



The solution to Eq. (J36J) is 



(41) 



(42) 



The explicit forms for P* tXX , P* >y y> and P\ xy can be found in the Appendix. The corresponding expressions for species 
2 are easily obtained by adequately changing the indices. This solution is still formal as we do not know the shear-rate 
dependence of A and the temperature ratios 71 and 72. These quantities must be consistently determined from the 
requirements 



£i7i 



p* 



(d - i)pr 



(43) 



£272 



2,xx 



(d-i)P 



2,yy 



72 



d 

1 - X171 



(44) 
(45) 



which follows from Eq. ([8]). Since the collision frequencies Lo rs are nonlinear functions of the temperature ratios, then 
it is not possible to get a closed equation for A or j r . Thus, one has to numerically solve the set of nonlinear equations 
(HP and (HU). 

Nevertheless, there are some limiting cases for which the problem can be solved analytically. For instance, in the 
case of mechanically equivalent particles (mi = 7712, <J\ = ^2,(^11 — «22 = 0:12), one gets that 71 = 72 = 1 and the 
partial pressure tensors P*+j can be written as 



A, 



x\ 



2,yy 

X2 



1 



1 + 2A' 



(46) 



p* 

l,xx 
Xl 



Pl xx 1 + 2dA 



■1-2 



1 + 2A ' 



(47) 



p* 
1,3:1/ 

Xl 



P 



X2 



(1 + 2A) : 



(48) 



where 



(49) 
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and A is the real root of the equation 



A(1 + 2A) 2 



d ' 



namely, 



A(a) = - sinh 2 
v ; 3 



- cosh" 1 [ 1 
6 



27. 



(50) 



(51) 



The parameter A governing the long-time behavior of the temperature can be easily obtained from Eqs. (|3ip and (|40l 

as 



A 



2n* 

_>* _ z " p * 
d v 

d + 2 



2d 



where use has been made of the result C* = (1 — a 2 )/ 2d. Equations ((46 
monocomponent gas @, Moreover, in the absence of shear field (a* 
viscosity function rf = rj* + r/^ where 



0), P* 



(B* n - C)(B* 2 - C) - B* 12 B* 21 
The temperature ratio 7 = 71/72 is determined from the condition 

xi _ Pl xx + {d- l)Pj vv 



(52) 

52]) are the same as those obtained for a 
= P* yy = x r ~/ r and the shear 

(53) 
(54) 



These results are consistent with those obtained for IMM in the Navier-Stokes regime [2l| . 

It must be remarked that, although the scaled pressure tensors P^ij reach stationary values, the gas is not in general 
in a steady state since the temperature changes in time. Actually, according to Eqs. (|3"Tj) and (fM)) . one gets 



-c 



2a* 



-V ' 



(55) 



Equation ([51))) shows that T(t) cither grows or decays exponentially. The first situation occurs if 2a* 2 rf > d(* . In 
that case, the imposed shear rate is sufficiently large (or the inelasticity is sufficiently low) as to make the viscous 
heating effect dominate over the inelastic cooling. The opposite happens if d(* > 2a* 2 rf . A perfect balance between 
both effects takes place when 2a* 2 rf = d(*. 



The expressions for the rhcological functions 77* (a*) and ^(a*) depend on many parameters: 
{xx, mi/m,2, 0"i/<72, a Ui a 22, Q!12) o*}- Obviously, this complexity exists in the elastic case as well [231 ] . so that 
the primary new feature is the dependence of 77* (a*) and x I'i(a*) on dissipation, on which we shall concentrate. Also, 
for simplicity, we take the simplest case of common coefficient of restitution (an = 022 = Q!i2 = &)• This reduces 
the parameter space to five quantities: {x\, mi/iri2, 01 /02, a , a*}. Before considering the rheological functions rj*(a*) 
and ^>i(a*), it is interesting to analyze the dependence of the temperature ratio T1/T2 on the shear rate. This 
quantity measures the lack of equipartition of the kinetic energy. Obviously, T\ = T2 for any value of the shear rate 
and/or the coefficient of restitution in the case of mechanically equivalent particles. Figure [1] shows the temperature 
ratio as a function of the (reduced) shear rate a* in a situation where the grains have the same mass density 
[(<Ji/(J2) d = mi/7712]. As is often the case in driven binary granular gases, the more massive particles have a larger 
granular temperature for moderate shear rates, while the reverse conclusion holds at high shears. The inelasticity 
parameter a and shear rate a* are here considered as independent, which in general results in an unsteady situation 
for the system. The intersection of a curve 71 (a*) for a given a with the steady state line shown by the continuous 
thick line, provides the shear rate a* corresponding to an exact balance between viscous heating and inelastic 
dissipation. For a* > a*, the temperature diverges, while it decays to in the opposite case. We note that even in 
the clastic case, the temperature ratio differs from unity, as a signature of non equilibrium behaviour. Only when 
the shear rate does vanish do we recover the equilibrium equipartition result (71 = 1). Fixing dissipation at a = 0.7, 
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0,8 - 

i I i I i I i — 

0,5 1 1,5 2 

a* 



FIG. 1: Shear-rate dependence of the temperature ratio 71 = Ti/T in three dimensions (d = 3), for an equimolar mixture 
(xi — 0.5), with cti/<T2 = 2, and mi/m2 = 8. Three values of the (common) coefficient of restitution are displayed, together 
with the stationary curve, also shown in Figure [7] (see text for details). 




FIG. 2: Temperature ratio 71 = Ti/T as a function of the reduced shear rate a* in three dimensions (d — 3). Here, a denotes 
the size ratio (ri/02 while ^ = mi/mj. 

Fig. [5] complements the results of Fig. [T]by showing the influence of mixture composition x±. The same qualitative 
trends are observed as in Fig. [TJ see the three upper curves. However, Fig. [5] also shows that changing the mass ratio 
(other parameters being fixed) can alter the results and lead to an increasing ratio T\/T with increasing a*. 



It appears that the non-linear shear viscosity and viscometric function exhibit a more robust behaviour with shear 
rate. It can be seen in Fig. [3] that those functions decrease with increasing a*. In this figure, the steady state values 
are shown with the continuous thick curve, and the maximum possible stationary shear rate is indicated by the squares 
(which correspond to a = 0). As above, for a given inelasticity, the intersection of the 77* (resp ^*) curve with its 
steady-state counterpart determines the steady-state value of the shear rate (resp normal stress difference). As in 
Fig. [1] three coefficients of restitution have been chosen in Fig. 02 and correspond to strongly inelastic (a = 0.5), 



10 



1 

- \ 




i_ 
\ 


1 






\ 

\ 


IX— U.J 




- 


\ 

- V 

\ 


■-■ ot=0.7 
■■■■ a=l 




- 


\ 

\ 


— steady state 






\ \ 
- \ \ 








\ \ 








\\ 































•V o 



0,5 



1,5 



FIG. 3: Shear-rate dependence of the reduced non-linear shear viscosity rj* for d — 3, an equimolar mixture (xi — 0.5), 
o"i/o"2 = 2, and mi/m.2 = 8, for three values of the (common) coefhcient of restitution a (same situation as in Fig. [T}. The 
inset shows "J* versus the reduced shear rate a* . The continuous thick curves are the loci of steady-state results (also shown in 
Fig. when dissipation is scanned in the admissible range a £ [0, 1]. The squares locate the terminal shear rate at maximum 
dissipation. 




FIG. 4: Non-linear shear viscosity rf and normal stress difference 'J* versus a*. The parameters are the same as for Fig. [2] 
and the different curve styles have the same meaning as in Fig. [2] 



moderately inelastic (a = 0.7), and elastic systems (a = 1). For completeness, we also show in Fig. @]the rheological 
functions corresponding to the parameter set of Fig. [2] The same qualitative trend is observed as in Fig. [3l We note 
in Figures [21 [3] and 0] the systematic trend that at high shear rates, the shear viscosity and normal stress difference 
become practically insensitive to the parameters specifying the state of the system, in particular dissipation. This 
observation is reminiscent of the single species phenomenology 0]. We also observe that the dependence on mixture 
composition x\ and mass ratio is more subtle (see e.g. the <J\j<yi = 2 and mi/7712 = 8 curves showing that the largest 
shear stress and normal stress difference occur in the equimolar case). Likewise, it is observed that for x\ = 0.5, 
ci/<72 = 2, the smallest shear stress and normal stress difference correspond to like masses (mi = 7712). We do not 
dwell on those effects since they are already present in the vanishing shear rate limit. A non-trivial effect of shear 
rate, however, is illustrated in Fig. [5j whereas at small a*, the viscosity function decreases with increasing a (at least 



11 



in the physically relevant range a > 0.7), increasing the shear rate leads to the opposite effect (see the inset of Fig. 
[5]). Enhanced dissipation then leads to smaller shear stresses. In Fig. [5j we have displayed the full possible range 
< a < 1 to show that even at small shear rates, an extreme dissipation can lead to a decreasing shear viscosity. We 
finally note here that those effects are absent for the normal stress difference, that appears to be a decreasing function 
of a, see Fig. [6] 





B. Model B 



In Model B the collision frequency vq(T) is an increasing function of temperature, and so the reduced shear rate 
a* is not constant. In order to have r/*(a*) and ^(a*) in Model B, one has to solve numerically the non- linear set 
PS)) , discard the kinetic stage of the evolution, and eliminate time in favor of a*(t) 0, Iu addition, it should 
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be remembered that the above functions depend on the temperature ratio, that is itself time dependent through its 
dependence on a* (t) . The above task in the case of a mixture is therefore a significantly more complex problem than 
in the monodisperse case. However, the results derived in the single gas case Q indicate that the influence of the 
temperature dependence of Uq on the rheological properties is quite small. We then restrict here our discussion to 
the steady-state solution for Model B. In this case, it is easy to see that the results obtained in the steady simple 
shear flow state are universal in the sense that they apply both for Model A and Model B, regardless of the specific 
dependence of v§ on T. 

Figure [7] shows the corresponding rheological functions and temperature ratio obtained at long times for a selected 
parameter set. The Figure also illustrates that in the steady state, dissipation and reduced shear are coupled (see 
the graph on the right hand side): to every value of a* is associated a given inelasticity, so that a* vanishes in the 
elastic limit a — * 1. This explains why the stationary temperature ratio converges to 1 in the small shear limit. For 
a given set of parameters 7711/7712, cr\/<J2, and mole fraction x\, the curves displayed are obtained by scanning all 
possible inelasticity parameters a from 1 corresponding to a vanishing stationary shear rate, to a = 0, which yields 
the maximum possible value of a* (e.g. 0.24 for x = 0.1 for the parameters used in Fig. [3 as can be seen in both left 
and right hand side graphs). The equimolar case results have already been shown as the continuous curves in Figs. [T] 
and [3J Interestingly, it can be seen that the normal stress difference can become an increasing function of a* (see the 
inset of Fig. [7]), whereas it is decreasing within Model A, when a is fixed and a* is changed. This is an illustration of 
the conflicting effects at work at the stationary point, when simultaneously increasing a* and dissipation (indeed, 'J'J 
decreases when a is fixed and a* increases while it increases when a* is fixed and a decreases). As far as the shear 
viscosity is concerned, the effects at work always lead to a decreasing function of a*, as with fixed dissipation within 
Model A. 




FIG. 7: (left) Plots of the temperature ratio 71 = T1/T2 and the non-linear shear viscosity 77* as functions of the (reduced) 
shear rate a* in the steady state for d — 3, 0-1/02 = 2, mi/7712 = 8, that corresponds to grains of same mass density. Three 
different compositions are displayed. The inset shows the normal stress difference ip* = rp* versus reduced shear rate. The 
right-hand side graph shows the a-dependence of the (reduced) shear rate a* for the above three systems. 



V. CONCLUSIONS 



We have determined the rheological properties (shear stress and normal stress difference) for a binary granular 
mixture in a uniform shear flow. The problem has been addressed in the framework of the Boltzmann equation 
with Maxwell kernel. An important advantage of such a model -compared to the more realistic inelastic hard sphere 
kernel- is that the collisional moments can be obtained exactly, and do not require the explicit knowledge of the 
velocity distribution function of both species (only low order moments are required). It is important here to stress 
that in spite of its approximate nature, IMM have been shown to fare very favorably against Monte Carlo simulations 
of inelastic hard sphere mixtures [H, [22| . We therefore expect the trends reported here to be realistic and of practical 
interest. 

In the case of IMM, there is a characteristic frequency v$ that can be chosen freely. To study the rheological 
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properties, two classes of models have been introduced. In Model A, v$ is taken independent from the total kinetic 
temperature T, while in Model B, v scales like T q , with q = 1/2 to reproduce the hard-sphere behaviour. In this 
respect, Model A fundamentally differs from models with q ^ in that it allows to disentangle the effects of dissipation 
(embodied in the coefficients of restitution a rs ) from those of imposed shear (embodied in the reduced shear rate a* 
defined from the actual shear rate a by a* = a/vo). Indeed, whenever q 0, the system reaches at long times a steady 
state where viscous heating compensates for collisional cooling. In this steady state, a key point is that the reduced 
shear rate is enslaved to the inelasticity, so that the problem is inherently non-Newtonian [1 01 ] (in other words, it is not 
possible to decrease the reduced shear rate a* by decreasing a, since then collisional dissipation will be more efficient 
and lead to a smaller granular temperature, so that a* is finally unaffected). On the other hand, within Model A, 
inelastic dissipation and viscous heating generically do not compensate, so that the temperature of the system either 
grows without bounds, or decreases to 0. If, however, the reduced shear rate is adjusted to the precise value that 
is reached within Model B for a given parameter set, Model A admits a steady state (which is then identical to its 
Model B counterpart). 

Outside the particular steady state point and within Model A, the scaled pressure tensor reaches well-defined 
stationary values which are non- linear functions of a* and the coefficients of restitution a rs ■ This allows for a clear- 
cut definition of the (reduced) non-linear shear viscosity rf and normal stress difference ^\ (within Maxwell models, 
only one stress difference is non- vanishing, whereas there are two such quantities for the hard sphere kernel). These 
quantities have been computed for various parameters characterizing the mixture, and it has been found that when 
a* is large enough (say a* > 1), rj* and "F^ are insensitive to dissipation, mixture composition, mass ratio, and size 
ratio, but only depend on the shear rate. Both quantities decrease when a* increases. At smaller shear amplitudes, 
the detailed parameters of the mixture becomes relevant and in particular we have found -for physically relevant 
dissipation parameters- that the shear viscosity and normal stress difference increase when dissipation is increased. 
The dependence on mixture composition and mass ratio appear more subtle, and are non-monotonous. Of particular 
interest also is the temperature ratio, which has been seen to depend on the parameters of the problem in a complex 
fashion. 

Future developments include the study of the tracer limit, segregation of an intruder by thermal diffusion, together 
with a more general derivation of generalized transport coefficients from a Chapman-Enskog-like expansion [25[ . Work 
along these lines is in progress. 



Acknowledgments 



This work was initiated during a visit of V.G. to the Laboratoire de Physique Theorique et Modeles Statistiques, 
Universite Paris-Sud. He is grateful to this institution for its hospitality and support. The research of V.G. has been 
supported by the Ministerio dc Education y Ciencia (Spain) through grant No. FIS2007-60977, partially financed by 
FEDER funds and by the Junta de Extremadura (Spain) through Grant No. GRU08069. 



APPENDIX A: PARTIAL PRESSURE TENSORS 

The expression of the relevant elements of the partial pressure tensor are given by 

Px<*x = ^3 ( GA2s i2 + 2a* 2 G5 12 [3A 2 + B% + B* 12 B* 21 + B* 22 (B* 22 + 3A) + B* n (B* 22 + 3A)] 

-FA 2 (B* 22 + A) - 2a* 2 F [B* n B* 12 B* 21 + B* 12 B* 21 (2B* 22 + 3A) + (B* 22 + A) 3 ] } , (Al) 

PT vv = GBl2 ~ F } B22+X \ (A2) 



where 



Pi.xy = i F [ B ^B 2 i + (B 22 + A) 2 ] - GB 12 (B n + B 22 + 2A)} , (A3) 
A = B* 12 B* 21 - (B* u + X)(B; 2 + A), (A4) 



F = A* n pl + A 12P * 2 , 



(A5) 



14 



G = A* 2 p* 2 +A nP * v (A6) 

In the above equations, A* s — A rs /v and B* s = B rs /v where A rs and B rs are given by Eqs. P7|) - ([5Q)) . The 
expressions for P| i? - can easily obtained from (|Al|) - (|A3l) by the adequate changes of indices. 
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